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Abstract 

We examine the Pioneer anomaly - a reported anomalous acceleration affecting the 
Pioneer 10/11, Galileo and Ulysses spacecrafts - in the context of a braneworld scenario. 
We show that effects due to the radion field cannot account for the anomaly, but that a 
scalar field with an appropriate potential is able to explain the phenomena. Implications 
and features of our solution are analyzed. 

1 Introduction 

Studies of radiometric data from the Pioneer 10/11, Galileo and Ulysses have revealed the 
existence of an anomalous acceleration on all four spacecrafts, inbound to the Sun and with a 
(constant) magnitude of a a — (8.5 ±1.3) x 10 -10 ms~ 2 . Extensive attempts to explain this 
phenomena as a result of poor accounting of thermal and mechanical effects, as well as errors in 
the tracking algorithms used, have shown to be unsuccessful [1] , despite a recent claim otherwise 

[2]- 

The two Pioneer spacecrafts follow approximate opposite hyperbolic trajectories away from 
the Solar System, while Galileo and Ulysses describe closed orbits. This, together with the 
fact that the three designs are geometrically distinct, explains the lack of an "engineering" 
solution for the anomaly. However, it prompts for a much more intriguing question: what is 
the fundamental, and possibly new, physics behind this anomaly? 

To answer this, many proposals have been advanced. The range of ideas is quite diverse and 
we mention some of them: Yukawa-like or higher order corrections to the Newtonian potential 
[3]; a scalar-tensor extension to the standard gravitational model [4]; Newtonian gravity as a 
long wavelength excitation of a scalar condensate inducing electroweak symmetry breaking [5] ; 
interaction of the spacecrafts with a long-range scalar field coupled to gravity [6, 7] ; an inverse 
time dependence for the gravitational constant G [8]; a length or momentum scale-dependent 
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cosmological term in the gravitational action functional [9]; a 5- dimensional cosmological model 
with a variable extra dimensional scale-factor and a static external space [10]; a local Solar 
System curvature for light-like geodesies arising from the cosmological expansion [11]; similarly, 
a recent work argues that the reported anomaly is related with the cosmological constant at the 
scale of the Solar System [12], even though this would lead to a repulsive force; an interaction 
with "mirror gas" or "mirror dust" in the Solar System [13]; a superstrong interaction of 
photons or massive bodies with the graviton background, yielding a constant acceleration, 
proportional to the Hubble constant [14]; expansion of solid materials on board deep space 
probes and contraction on Earth due to a curved stress field arising from repetitive tidal action 
[15]; an expanded PPN-framework so to incorporate a direct effect on local scales due to the 
cosmic space-time expansion [16]; a result of flavor oscillations of neutrinos in the Brans-Dicke 
theory of gravity [17]; a theory of conformal gravity with dynamical mass generation, including 
the Higgs scalar (capable of reproducing the standard gravitational dynamics and tests within 
the Solar System, and yet leaving room for a Pioneer-like anomaly on small bodies) [18]; a 
gravitational frequency shift of the radio signals proportional to the distance to the spacecrafts 
and the density of dust in the intermediate medium [19]; resistance of the spacecrafts antennae 
as they traverse interplanetary dust [20]; a gravitational acceleration a oc r -2 for a constant 
a a = 10~ 10 ms~ 2 and a oc r" 1 for a C ao [21]; clustering of dark matter in the form of a 
spherical halo of a degenerate gas of heavy neutrinos around the Sun [22] - amongst a few other 
suggestions put forward. It is interesting to mention that in higher-curvature theories of gravity 
where the gravitational coupling is asymptotically free and which have been much discussed in 
the context of the dark matter problem [23, 24, 25], a stronger gravitational coupling is expected 
on large scales and hence, at least in principle, to a Pioneer-like anomalous acceleration. 

In this work we consider the Pioneer anomaly in the context of the braneworld scenario. 
We use the Randall-Sundrum model and variations to show that gravitational effects such as 
the one due to the radion field cannot explain the anomaly. We argue that a scalar field with 
a suitable potential implies that geodesies in this theory exhibit an extra constant attractive 
acceleration. 

2 Braneworld Theories 

A quite new range of possible scenarios arise in the context of braneworld theories. In these, 
one assumes our Universe to be a 3-dimensional world-sheet embedded in a higher dimensional 
bulk space. Considerations on the symmetries of the brane and its topological properties are 
then taken to constrain the evolution of matter on the brane and gravity on the brane and in 
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the bulk. 



Braneworld theories are a fast developing trend in cosmology [26, 27, 28, 29, 30], whose 
main feature is to allow for a solution for the hierarchy problem, whether the typical mass scale 
of the bulk is comparable with the electroweak breaking scale, M EW ~ TeV. 

In this work, we shall consider the Randall-Sundrum braneworld model and some variations 
[26, 27]. One admits a scenario with two 3-branes embedded in (4+l)-dimensional space: a 
positive tension brane situated at z = and a negative tension one at z — z c . The Ansatz for 
the metric takes the form 



which is a solution of the 5-dimensional Einstein's equations and preserves Poincare invariance 
on each brane. The constant k is a fundamental quantity of the theory and typically takes 
values of order k ~ Mpi, the 4-dimensional Planck mass, which is dynamically generated from 
the bulk space Planck mass M 5 , 



This relation is obtained from the derivation of a 4-dimensional effective action by integrating 
the fifth dimension away. Notice that Mpi depends weakly on the second membrane position, 
in the large kz c limit. Physical masses on the positive brane, however, scale with this distance 
through the relationship m = e~ kZc mo (a solution for the hierarchy problem is obtained for 
kz c ~ 15). In traditional compactification schemes, the "warp" factor e~ 2kz in the metric is 
absent, and hence integration over the fifth dimension yields only the volume of the bulk space, 
V n ; for a n-dimensional compact space one obtains M Pl = M n+2 V n . 

3 Braneworld Scenarios for the Pioneer anomaly 

As a first attempt to explain the Pioneer anomaly within the context of braneworld theories, 
one could resort to the appearance of a tower of Kaluza-Klein massive tensorial perturbations 
to the metric. Three problems arise: all gravitons are ordered according to their mass, so 
that one cannot freely specify the range of one of them without affecting the whole tower. 
Most braneworld models consider one first light mode with cosmological range, and all ensuing 
modes to have sub-millimeter range. It is difficult to introduce an intermediate scale without 
abdicating from one of the two desired extreme cases. 

The second problem refers to the fact that any Yukawa gravitational potential would affect 
all bodies within range, independently of their mass (as expected from the Equivalence Princi- 



d s 2 = e~ 2kz g ilv dx il dx v + dz 2 , 



(1) 




(2) 
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pie). This is in direct contradiction with the lack of an observed "anomalous" acceleration for 
the planets in the Solar System. 

Thirdly, it has been shown that an Yukawa potential fitted to the observed effect would 
reveal the presence of a graviton with range L ~ 200 AU and a negative coupling a ~ — 10~ 3 
[3]. One can adjust the exact values, since they belong to a solution curve L = L(a) (the one 
presented is a rather "natural" choice, since |a| <C 1 and L ~ 100 AU). However, a is always 
negative. Since the coupling of modes is proportional to the normalization factor of their wave 
functions in the Kaluza-Klein reduction scheme, this would imply in a graviton with negative 
norm - which is unattainable within current braneworld theories, leading to instabilities at a 
quantum level. This latter issue shows that quantized tensor excitations are not the cause of 
the Pioneer anomaly. 

If one considers instead a single membrane, then the periodic boundary condition disap- 
pears and the modes are no longer quantized. In the standard Randall-Sundrum scenario [27], 
this amounts to a modification of gravity at large distances, with the gravitational potential 
presenting a 1/r 2 behaviour, typical of the five-dimensional propagation of gravity in the bulk. 
This behaviour also appears in models where the modes are quantized and dense, and their 
mass spectrum approaches a continuous distribution [30]. 

A more elaborate model was suggested in Ref. [29], which shows a similar behaviour of the 
gravitational potential. At large distances it also goes as 1/r 2 , while at intermediate distances 
it contains a repulsive logarithmic term. In both cases (and more exotic ones, see Refs. [30, 31] 
and references therein), a continuous spectrum of Kaluza-Klein excitations can be shown to be 
unable of explaining the Pioneer effect. 

An alternative solution could be that the Pioneer anomaly reflects the influence of the radion 
field f(x), a scalar perturbation of the metric corresponding to translational zero modes, related 
to relative motion of the two branes. Since its moduli is usually stabilized due to some ad hoc 
potential whose effect is superimposed on the usual warped metric (see Ref. [32] and references 
therein), one could view the radion as an additional field on each brane. We discuss its effect 
on a test particle and verify to which extent it could provide an explanation for the Pioneer 
anomaly. 

Following Ref. [33], the induced metric on the positive (z = 0) brane which includes this 
perturbation is, in Gaussian coordinates (g zz = — 1 and g z = 0), 

h; u = % „[l-2kf(x)] + ±-f^ , (3) 
where f(x) is constrained by the requirement that in the vacuum 
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and 



□ 2 / = o 



(4) 



□ 2 = (v X °U) a = V Xu fM ~ V Xa TLU (5) 

is the 4-dimensional d'Alembertian. This arises from the Israel junction conditions imposed 
on each brane so to ensure the Z & symmetry. The metric i] Xa is the unperturbed one on the 
positive brane. In the presence of matter, Eq. (4) is not homogeneous and is related to the 
trace of the energy-momentum tensor of the matter distribution. 

We now search for a solution /(r) that is both static and spherically symmetric. We work 
in spherical coordinates where /' denotes f\ r , and f t g = = f\ t = 0. In the weak field limit, 
one has rj o = 1 + 2U and r\ rr = — (1 — 2U), where U = —GM Q /r = —C/r is the Newtonian 
gravitational potential; in units where c — h — 1, C — 1.5 km. Thus, one obtains 

U 2 f = -f" --f' = , (6) 

r 

neglecting the term due to the curvature, which is proportional to U'(r)f'(r). 

A simple solution is given by f(r) = k~ 1 (A/r + B), with 5 C 1 a dimensionless constant, 
[A] = L and A/r C 1 within the desired range. Hence, the induced metric becomes 



h$ = [l + 2U(r)][l-2kf(r)] 
hX = -[l-2U][l-2kf(r)]+k- 1 f"(r) . 



(7) 



We next proceed by considering the spacecrafts as point particles. Thus, they follow time- 
like geodesies of the obtained vacuum metric. The acceleration is then 



-r 00 = h rr d r h 00 = -^[i + 2U] 



1 - 2kf - k~ l f"(r) 



d r [(l + 2U)(l-2kf)} 



~ -[1 + 2U] [1 + 2k f + k- l f"]d r [U-kf- 2kUf)} 
C + A -2A 2 + 4AC + 2C 2 4A(3A - C)C 



(8) 



2A(A + C + 8AC 2 k 2 ) 4AC(3A + C) 16A 2 C 2 



k 2 r J 



+ 



k 2 r 6 



k 2 r 7 



and one can see that no constant term arises, rendering this approach unsuitable to account 
for the discussed anomalous acceleration. 
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Another possible explanation for the Pioneer anomaly could be related with a bimetric 
theory exhibiting Lorentz symmetry breaking. A proposal along these lines was discussed in 
Ref. [34] and is briefly presented in the Appendix A. 

In principle, a bimetric theory could arise from induced effects on the 3-brane of gravity in 
the bulk. However, as shown in Ref. [35], Goldstone modes resulting from the spontaneous 
breaking of coordinate diffeomorphisms carry no extra degrees of freedom due to the positions 
of the branes in curved spacetime, but rather manifest themselves in the form of an extra 
field in the energy-momentum tensor, involving the induced metric 7 mn = d m x IJ 'd n x u g l j il/ (x). 
Thus, although invariance under general coordinate transformations does not hold along the 
fifth dimension, it is a symmetry on each brane. Since the resulting singular energy-momentum 
tensor cannot be treated as a "mass" term for the graviton (at least at the linear level), massive 
gravitons do not arise due to the interaction of the brane with gravity in the bulk. It is argued 
that this is crucial to preserve the r~ 2 long-range behaviour of gravity. 

Nevertheless, it should be pointed out that long-range modifications to Newton's law have 
been the focus of many braneworld related proposals aiming to explain the accelerated ex- 
pansion of the Universe (see e.g. Ref. [28] and references therein). It is important to realize 
that the above arguments are valid in the absence of topological obstructions on the brane. 
Hence, no extra degrees of freedom arise when the metric is taken as a dynamic variable un- 
der local, "small", diffeomorphisms. However, under large gauge transformations [36], where 
the gauge parameters assume different values in non-connected asymptotic regions and have 
different global behaviour, solitonic-type solutions are admitted. 

These solutions are deformations, h^ v , of a particular metric solution gfy of the Einstein's 
equations with a source term given by the extra field arising from the breaking of diffeomorphism 
invariance along the fifth direction, such that the deformed metric g$ + is a solution of 
the same equation. We have verified that known solitonic-type solutions do not give rise to a 
constant acceleration term [37, 38] that could explain the Pioneer anomaly. 

In what follows, we shall derive the general behaviour of a scalar field with a potential and 
calculate its effect on the motion of a test particle. 

4 Scalar Field Coupled to Gravity 

As a possible explanation for the Pioneer anomaly, we consider the effect induced by the 
presence of a scalar field with dynamics driven by a potential V(<f>) oc — _Q (r), with a > 0. 
The form of this potential closely resembles that of some supergravity inspired quintessence 
models [39, 40, 41]. In the context of the braneworld scenario, the quintessence potential has 
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the form V oc 4>~ a (t), with 2 < a < 6 [42]. Notice that, in this proposal, there is a spatial and 
not a time dependence. Also, the sign of our potential is reversed, so to yield a static attractive 
acceleration. 

As usual, we assume a small perturbation to the Minkowsky metric and solve it in terms of 
the energy-momentum tensor of the field. 
The metric can be written as 

fiV = + V ( 9 ) 
where in the weak field limit (in spherical coordinates) 

(rfiiw = diag(l + 2U(r),-l +2U(r),-r 2 ,-r 2 sin 2 6) , (10) 

and 

V = dtag(f(r), -h(r), -h(r)r 2 , -h{r)r 2 sin 2 6) . (11) 

Notice that the bimetric character arises from the assumption that the field <fi expresses the 
effect of the induced metric arising from the spontaneously broken diffeomorphisms in the 
curved spacetime. 

The Lagrangian density of the static scalar field takes the form 

^ = lvMd»d> - v(4) = \ nrr (0f + A 2 r a , (12) 

where A is a constant. The scalar field obeys the equation of motion 

rf*+T=°- (13) 

which yields, neglecting a term proportional to U'(r)(f) f (r) 

4>"(t) + -4>'{r) = aA 2 ^- 1 , (14) 
r 

and admits the solution 

2 

0(r) = ((2 + a^^L-A^J ^ = . (15) 

Notice that this solution is singular at r = 0. Its regularization is discussed in Appendix B. 
Thus, in terms of r, the potential and gradient terms are given by 



V((f)(r)) = -A 2 (3 a r~^ , 
7 



(16) 



and 



W'll'^I^A---^)^*)] ■ (17) 



1 / 2 .o ( a \ _ _^o L / a 

2 

The Lagrangian density is given in the Newtonian limit by 



£ * = -4T^ l ' Wr)) = 4T^^* ■ (18 > 

The energy-momentum tensor of the scalar field is obtained by the expression 



Tp, = d^d^ - 7)^X4, , (19) 
so that its components are given by 



4A 2 



Too — 




T — 

rr 


0'(rf - 


Tee = 




T — 
j- pip 





'4 + a 
, 2A 2 



• * v '4 + a 

2 4A 2 aa 

4 = r- (3 r 2 +« 



4 + a 

4A 2 



„2 „• 2/ 



(20) 



4 + a' 

where we have assumed that the spatial perturbation to the metric is very small, h(r) <C 1 



The trace is given by 



2A 2 

T = T a =r] ^ T ^ = _ {8 + a) __ (3 a r -^ (21) 

We now turn to the linearized Einstein's equations: 

^VV = «(^-^V r ) > ( 22 ) 
where k = 87rG. The — component is 

f»(r) + ^/'( r ) = 2« (l - ^) ^/Tr"^ (23) 
whose solution is, for q^2, 



and, for a = 2, 
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/(r) = ^_^ CKi09 (_L.) . ( 25 ) 

Notice that we have dropped the homogeneous solution of Eq. (23) of the form 1/r since it can 
be absorbed by the Newtonian term. Therefore, we obtain for the acceleration 

a r = -- 2 + (2 + a)A"^ a r'^(^--^—) . (26) 
r z \ 2 6 + ay 

For a = 2, it reads 



a r = -^-\l«— + \l«— ■ (27) 



C 3Ak 3ACk 
7 1 ~ \2~ + V 2~V 

Eqs. (26) and (27) indicate that a constant anomalous acceleration exists only for a — 2. 
The first term is the Newtonian contribution, and one can identify the constant term with the 
anomalous acceleration a A = 8.5 x 1CT 10 ms~ 2 , setting 

A = 4a/^— = 4.7 x 10 42 m" 3 . (28) 

The remaining term, proportional to r _1 , is much smaller than the anomalous acceleration 
for AC /r <C 1, that is, for r 3> 6 km; it is also much smaller than the Newtonian acceleration 
for r <C 2.9 x 10 22 /cm 100 Mpc, clearly covering the desired range. It can be shown that 
higher order terms affect the third term in Eq. (27) and yield a negligible correction to the first 
term. 

For consistency, we now show that h(r) 1. For that, we use the r — r component of Eq. 
(22) with a = 2: 

9 Ak 

h"(r) + Z -h'(r) = , (29) 

r V6r 

whose solution is h(r) = —Anr/2\/Q. This term is negligible for r <^ 5000 Mpc, confirming the 
validity of the approximation. One can improve the previous result of Eq. (27) to first order 
in h(r): 



(i) a A „ C 5 ^ACk fz An / /2 \ (Ak\ 2 

a (1) = f— ~ a A (l-/i(r)) = — 7 + \/ 1 - \\-ACk - — r , (30) 

l + /i(r) V v;; r 2 6V2r V24\ v V3 j Ui 1 ; 

which, asides from the small correction to a A and a 5/6 factor in the r~ l term, introduces a 
linear term which is negligible for r< 5 x 10 29 Mpc. 

Thus, we see that an anomalous acceleration is a clear prediction of our model within the 
Solar System. Hence, an hypothetical dedicated probe to confirm the Pioneer anomaly [43, 44] 
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does not need to venture into too deep space to detect such an anomalous acceleration, but 
just to a distance where it is measurable against the regular acceleration and the solar radiative 
pressure, actually approximately from Saturn onwards. 

It can be shown that the potential V(4>) oc — 0~ 2 is the only one that can account for the 
anomalous acceleration. Indeed, from Eq. (23) we can write in the same approximation 

^[r 2 f'(r)Y = XKV(r) m (r) , (31) 
where x is a dimensionless factor depending on the potential (x = —2 for a — 2). Hence, 

/» = ^ / r 2 V(r)r) Q0 (r) . (32) 
Since, from Eq. (27), the anomalous acceleration is given by 

or = -\f\r) = / V(ry mo (r) , (33) 

and r/oo = 1 — 2C/r, one concludes that the potential must have the form V(r) = V(<f>(r)) = 
—A/r or V(r) = B, where B is a constant. The constant potential solution also provides an 
anomalous acceleration, but yields a conflicting linear term: 

f (const) = -2BCnr + Bnr 2 /3 (34) 

and hence 

C Bur 

a r (const) = — o ~ BCk + ■ ( 35 ) 

r z 3 

By identifying the constant term with the anomalous acceleration, one immediately finds 
B = 0.03 kg m~ 3 ; as a result, the linear term dominates the constant acceleration for r 3> 
1.5 km. Hence, we consider the case x — ~ 2 or a = 2. 

Indeed, the equation of motion (13) can be written as 

V 2 = ^[rV(r)]'=^ > ( 36 ) 
r A (p'[r) 

and therefore 

0'(r)[rV(r)]' = A . (37) 

Thus, the only real solution of this differential equation is given, up to a constant factor, by 
y/r. Hence, r oc 2 (r) and V(<f>) oc — _2 (r), as argued. 
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With a — 2, the energy of this scalar field grows logarithmically. This can be seen as a 
manifestation of a global symmetry breaking, as in the case of global cosmic strings, where the 
same behaviour is found [46]. 

Notice that these results apply only to the case of bodies with negligible mass, such as the 
considered spacecrafts. In the case of planets, endowed in general with rotation and spin, their 
dynamics are described by the solution of Eq. (22) with the appropriate energy-momentum 
tensor T^ u = T^ u + T^ net , with T^ v <C T^ net . Hence, no anomalous acceleration is expected 
for celestial objects. 

Furthermore, in what concerns the nature of the source of the scalar field, we assume it is 
the Sun and that matter couples with the scalar field in the following way: 

C-int = <PY1 fSi^i » ( 38 ) 
i 

where, ipi stands for different fermion species and /j are couplings constants. 

We now look at a possible breaking of the Equivalence Principle. It can be shown (see 
[47, 48] and references therein) that changes to geodesic motion arise (to first order) from 
spatial variations of mass, proportionally to m'(r)/m(r). Hence, one must look at changes 
induced on the mass of a test particle. These will occur if the scalar field acquires a vacuum 
expectation value, such as in the standard case of the Higgs mechanism, or in the cosmologically 
relevant cases ([47]). 

In our proposal, the potential for the scalar field is negative and monotonically increasing. 
Hence, the effective potential V e //(</>) = V(4>) + (pj^ifi^i, where n» is the density of different 
matter species coupled to the scalar field, does not develop a minimum. Therefore, no mass 
changes as well as no violations of the Equivalence Principle are expected. 

Finally, we look at the propagation time delay induced by this bimetric theory. If it is 
not negligible, then the value of the anomalous acceleration should be corrected, since the 
radiometric data which supports it is based on the Doppler effect [1]. 

For simplicity, we consider only the worst-case scenario: a test particle travelling at a 
constant velocity of v — 10~ 5 c (approximately the current velocity of Pioneer 10/11), following 
a linear path away from the Sun and at a distance r from it. If the light signal is emitted with 
a proper period T , its relation with the period recorded on Earth, T, is 



T = I y/[l + f(r(t))]-v 2 dt~ [ 

JTq JTq 



Vl-v 2 + 



f(r) 



dt 

v 



Ar 



y/l - v 2 + J -Ak 



i^-C log 



vT 



dr 

v 



(39) 
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where Ar = vT . This yields 




3 Akt 



) 



To , 



(40) 



where we have used that r ^> 1.5 km and r ^> Ar. Hence, the bimetric effects on the time 
delay can be safely disregarded for r <C 3000 Mpc. 



In this paper we have investigated the possibility of explaining the Pioneer anomaly within 
the framework of braneworld scenarios. We have eliminated both "new" tensorial (massive 
gravitons) as well as scalar (radion) degrees of freedom as candidates for a solution. We found 
that the anomalous acceleration could be due to the presence of a negative potential scalar 
field, with a potential V oc — <fi~ 2 (r) similar to some supergravity inspired quintessence models. 

Notice that the approach considered here, contrary to naive thinking, has no implications 
for the puzzle of the rotation curve of the galaxy. Indeed, assuming a galaxy to be virialized, 
one can describe the rotation of a layer at a distance r of the galactic core as v 2 (r) = GM(r)/r, 
where M(r) is the total mass inside the layer. Observation shows that v 2 (r) displays a steady 
rise until a threshold of about 10 kpc (~2 x 10 9 A£7), and a constant plateau from there (see 
eg. Ref. [45]). This leads one to model M Ga i(r) ~ r and hence to postulate the presence of 
dark matter. 

Starting from Eq. (27), we can derive a different expression for r 2 (r): 



where the superscript * refers to galactic values. This curve does not describe the observed 
data and would lead one to model M(r) as 



This clearly differs from usual dark matter models due to the higher order terms in Eq. 
(42). However, likewise the case of planets, the test masses here cannot be viewed as point-like 
objects, and hence have to be treated with its energy-momentum tensor. Therefore, we are 
lead to conclude that the origin of the flattening of the rotation curves of galaxies does not 
have its roots in the induced bimetric theory. 



5 Conclusions and Outlook 




(41) 




(42) 
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6 Appendices 



6.1 Appendix A. 



A bimetric theory of gravity was first proposed by Rosen [49]. Its action is given by 



S 



1 



/ 




(43) 



Q4ttG 



where the vertical line | denotes covariant derivation with respect to the background metric 
r\^ v only, and Cm is the matter Lagrangian density. The resulting equation for the dynamical 
gravitational field is given by: 



where T = T^ v g^ u and ^ is the d'Alembertian operator with respect to r]^ v . Notice that 
the momentum-energy tensor couples only to the dynamical metric g^ v . We can always 
choose coordinates in which (77^) = diag(— 1, 1, 1, 1), the Minkowsky metric, and (g^ u ) = 
diag(—Co, c±, c±, Ci), where Co and C\ are parameters that may vary on a Hubble H~ l timescale 



This theory explicitly breaks Lorentz invariance. This is better understood by resorting 
to the Parametrized Post-Newtonian (PPN) formalism: a systematic expansion of first-order 
1/c 2 terms in the Newtonian gravitational potential and related quantities [50]. It turns out 
that any metric theory of gravitation can be classified according to ten PPN parameters: 
7, (3, £, cki, «2, «3, Ci> C2, C35 C4- These are the linear coefficients of each possible first-order term 
(generated from rest mass, energy, pressure and velocity), and relate a particular theory with 
fundamental aspects of physics: conservation of linear and angular momentum, preferred-frame 
and preferred-location effects, nonlinearity and space-curvature per unit mass, etc. 

Einstein's General Relativity, the most successful theory up to date, exhibits a set of PPN 
parameters with (5 — 7 = 1, the remaining being equal to zero. Rosen's bimetric theory has 
(5 — 7 = 1, but a non- vanishing a 2 = co/ci — 1 coefficient. This indicates that the theory is semi- 
conservative (it lacks angular momentum conservation) and exhibits preferred-frame effects: 
Lorentz invariance is broken and the Strong Equivalence Principle does not hold. It is worth 
pointing out that the breaking of Lorentz invariance has been recently very much discussed. 
Indeed, possible signatures of the breaking of this symmetry arise from ultra-high energy cosmic 
rays with energies beyond the Greisen-Zatsepin-Kuzmin cut-off, Eqzk — 4 x 10 19 eV, (see 
Ref. [51] for a discussion on the astrophysical aspects of the problem), from the observation of 
gamma radiation from faraway sources with energies beyond 20 TeV, and from the longitudinal 
evolution of air showers created by ultra-high energy hadronic particles and which seem to imply 



2K(g/ V ) l ' 2 (T, u - l -g, v T) 



(44) 



[49]. 
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that pions are more stable than expected (for an update see [34, 52] and references therein). 
Of course, Lorentz invariance holds with great accuracy as observed deviations are quite small, 
5 < 3 x 1CT 22 [53] from direct measurements, and even smaller from the study of ultra-high 
energy cosmic rays S ~ 1.7 x 1CT 25 [54, 55]. 

By linearizing Eq. (44) in the vacuum, one obtains the wave equations for weak gravitational 
waves, whose solution is a wave propagating with speed c g = \Jc\/cq. Thus, a 2 measures 
the relative difference in speed (measured by an observer at rest in the Universe rest frame) 
between electromagnetic and gravitational waves, inducing a time delay in the propagation of 
light signals [50]. This effect has been claimed to be put under test in a recent experiment 
using the close celestial alignment of Jupiter and the quasar J0842 + 1835. This analysis yield 
Cg/c = 1.06 ± 0.21 [56], corresponding to a 2 = —0.11 ± 0.35. Note, however, that this result 
is still controversial, and an alternative interpretation suggests that it sets instead a limit on 7 
and cki, actually 7 = 1 and a.\ = 0, being unrelated to the velocity of propagation of gravity 
[57]. 

A rigorous study of deviations between the Sun's spin axis and the ecliptic has led to the 
experimental constraint \a 2 \ < 1.2 x 10~ 7 [58]. Improving this bound as well as finding new 
means of verifying its implications are clearly of key importance. Interestingly improvements on 
the measurement of Sun's oblateness (and ensuing spin) as well as on the PPN parameters /3, 7 
and the combination r\ = 2 — f3 + 27 are on the list of objectives of the ambitious BepiColombo 
mission to Mercury [59]. 

Due to their small mass, self-gravitation is also absent for the considered spacecrafts. Thus, 
these can be regarded as particles, which enables us to calculate their acceleration by simply 
computing the time-like geodesies of the full metric, h^ v = rj^ + g^ u . 

In the weak field limit, v <^ c, one has 

a 1 ~ -Pa, = \h iX d x h 0Q , (45) 
from which one can identify the radial anomalous component of the acceleration: 

a A = Cl W r U - V r c . (46) 

It can be immediately seen that, if Cq and C\ are homogeneous in space, the derived anoma- 
lous acceleration is not constant, which contradicts the observation. Therefore, we assume 
these parameters depend on the distance to the Sun, that is, c = c (r), c\ = Ci(r). Given the 
constraint 
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ot 2 = 



c (r) 
ci(r) 



< 4 x 10~ 7 , 



(47) 



we assume that both coefficients have the same r-dependence, so that a 2 is homogeneous. We 
consider the choice 



c = D r , Cl = Fr , (48) 

with D, F > so that the resulting anomalous acceleration is inbound. According to Ref. [58], 
\D/F — 1| < 4 x 10~ 7 and hence D ~ F. Note that D cannot be exactly equal to F, since this 
implies that a 2 = 0. 

Substituting Ansatz Eq. (48) into Eq. (46), we find 



D 

aA = ~^2 



D r 



D 
2~ 



2C „ 
1 Dr 

r 



(49) 



Hence, D = 2a A = 1.9 x 10 26 m 1 and we see that the distance-dependent contributions to 
a a are negligible for r lying in the interval 

[2C, D- 1 } = [3 km, 3.5 x 10 14 AU] , (50) 

which is consistent with the fundamental assumption that c , C\ <^ 1. 

Unfortunately, this elegant solution for the Pioneer anomaly cannot be taken seriously, given 
the behaviour of Rosen's theory in what concerns gravitational waves. Indeed, in Ref. [60] , it 
is argued that Rosen's bimetric theory is fundamentally flawed, since it predicts dipole gravi- 
tational radiation beyond the limits measured from the pulsar PSR 1913 + 16 [61]. Moreover, 
the solution proposed by Rosen of considering a combination of retarded and advanced gravita- 
tional waves [62] implies in contradiction with the observed quadrupole gravitational radiation 
from binary pulsars [61]. Thus, Rosen's bimetric theory cannot be considered a viable theory 
of gravity. 

6.2 Appendix B. 

Equation (14) is not satisfied at the origin, the center of the Sun, as solution of Eq.(15) is 
singular at r = 0. However, this solution can be regularized in the standard way by introducing 
a source term in the Lagrangian density (12). This can be performed by considering a thin 
shell model, dividing the space into two regions, r > r and r < r , with r a positive constant. 
Assuming that the Pioneer anomaly is verifiable everywhere in the Solar System, it follows that 
r is smaller than R Q , the Sun's radius - this assumption can be relaxed, since it is currently 
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not accessible to experiment, due to the suppression of the effect by the larger solar wind and 
Newtonian terms. 

Naturally, one assumes that the solution of Eq.(15), <p{r) = </>(r)+, is valid for r > r . For 
r < r , we consider the solution 

0(r) = ( j ) _(r) = C + ^A 2 C- a - 1 r 2 + C (r 4 ) , (51) 

where C is a regularization constant. Continuity of the scalar field at r = r implies that C 
is a solution of 4>-(r ) = + (r o ). Even though the full solution is now regular at r = 0, the 
derivative changes value at r = r as <f>' + (r ) ^ <j)'_ (r ) and, as before, Eq.(14) is not satisfied at 
r = r . A suitable solution requires adding to Eq. (14) the term (r ) — 4>'_ (r )) 5 (r — r ), 
which demands an additional term in the Lagrangian density Eq. (12): 

£^^ + 0(0' + (r o )-0'_(t o ))5(r-r o ) . (52) 

Hence, one can remove the "inner" solution 0_(r) by taking the r — * limit and keeping the 
source term in Eq.(52). 
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